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Abstract. The issue of the approximate isotropy and homo- 
geneity of the observable universe is one of the major topics 
in modern Cosmology: the common use of the Friedmann- 
Robertson- Walker [FWR] metric relies on these assumptions. 
Therefore, results conflicting with the "canonical" picture 
would be of the utmost importance. In a number of recent pa- 
pers it has been suggested that strong evidence of a fractal dis- 
tribution with dimension D ~ 2 exists in several samples, in- 
cluding Abell clusters [ACO] and galaxies from the ESO Slice 
Project redshift survey [ESP]. 

Here we report the results of an independent analysis of the 
radial density run, N(< R) °< R D , of the ESP and ACO data. 

For the ESP data the situation is such that the explored 
volume, albeit reasonably deep, is still influenced by the pres- 
ence of large structures. Moreover, the depth of the ESP survey 
(z 0.2) is such to cause noticeable effects according to dif- 
ferent choices of k-corrections, and this adds some additional 
uncertainty in the results. However, we find that for a variety of 
volume limited samples the dimensionality of the ESP sample 
is D w 3, and the value D = 2 is always excluded at the level 
of at least five (bootstrap) standard deviations. The only way 
in which we reproduce D » 2 is by both unphysically ignoring 
the galaxy k-correction and using Euclidean rather than FRW 
cosmological distances. 

In the cluster case the problems related to the choice of 
metrics and k-correction are much lessened, and we find that 
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ACO clusters have D ACO = 3 .07 ± 0. 1 8 and D ACO = 2.93 ±0. 15 
for richness class > 1 and ^. > 0, respectively. Therefore 
D = 2 is excluded with high significance also for the cluster 
data. 
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1. Introduction 

Recently Pietron ero and collaborat ors [h ereafter P&C] ( Pietro - 
nero et al. |1997[ Sylos Labini et al. |1996| , Baryshev et al. |1994| ) 
argued that a large number of samples give strong evidence 
that the distribution of clusters and galaxies is a simple fractal 
of dimension D ~ 2 [] on scales of several hundreds of Mpc, 
quite different from the "canonical" value of D = 3. 

This is an important claim because of its far reaching im- 
plications, and this issue, which dates bac k to th e beginning of 
the century (Charlier, 1908, see Peebles 1993 ), can be prop- 
erly addressed with long and careful analyses which, however, 
demand much deeper and better samples than those presently 
available in order to be able to explore an adequate range of 
scales (see e.g. Mc Cauley 1997, Hamburger et al. 1996). While 
2D constraints come from fluctuations in cosmic backgrounds 



It must be noticed that this value is different from the value D ~ 
1.6, sup ported by the same authors in the past years (Coleman & Pie- 



tronero 1992h, and D = 2 has been often discussed in the literature, 



because it would naturally arise locally from a geometrical do minan ce 
of planar structures, such as "pancakes", see f.i. Guzzo et al. (1991). 
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(see e.g. Peebles 1993), historically much work on this sub- 
ject has been done for decades in analyzing galaxy counts 
which, in a non evolving Euclidean Universe, should follow 
N(m) oc 10° 6m . Indeed, this behaviour has been observed in an 
intermediate magnitude range, m ss 15 -j- 17, (see e.g. Sandage 



1995). However, on the bright side one has to deal with magni- 



tude errors given from saturation of photographic plates and/or 
the relatively small volumes sampled which also reflect in un- 
certainties in locally the derived space galaxy density (cf Love- 
day et al. [1992[ Zucca et al. 1997, Maddox |l997| ), while on the 
faint side cosmological curvature and evolutionary effects be- 
come dominant and difficult to disentangle (Koo & Kron 1992, 
Ellis 1997). Therefore one really needs redshift information in 
order to avoid effects of time-space projections. 

In this paper we limit ourselves to an independent check 
on two of the samples discussed by P&C, the ESP galaxies 
(Vettolani et al. 119971) and the ACO clusters (Abell, Corwin & 



Olowin 1989), without touching upon the very general issue 



of a possible fractal distribution of the matter in the universe 
and its consequences: the interested reader can consult Peebles 



(1993), Coleman & Piet ronero (1992), Stoeger etal. (1987), 



Ehlers & Rindle r dl987p , Szal ay & Schramm (1985), Luo & 
Schram m ( |l992| ), Mc Cauley (\99l\ , Buchert ( jl997| ), Guzzo 
( 1997 ) and references therein. 



However, if the evidence were present in the data at the 
highly significant level claimed by P&C, the simple analysis 
presented here should be more than adequate in confirming the 
D = 2 claims. 

In section 1 we present the formalism, in section 2 the re- 
sults from ESP galaxies, in section 3 the results from the ACO 
clusters, and finally in section 4 the conclusions. 

2. Method 



In a previous discussion (Scaramella et al. |1991j) on the ori- 
gin of the Local Group velocity with respect to the Cosmic 
Microwave Background frame, it was pointed out the rele- 
vance of the depth's behaviour of clusters' "monopole" and 
"dipole" with respect to the issue of appr oach to homogeneity 
of the matter distribution (see Sylos Labini 1994 for contrasting 
views). 

The method, adopted also by P&C, is very simple: basi- 
cally, after selecting a volume limited sample of objects, one 
considers the quantity 



obtained by averaging on all points of the sample (for which 
several methods have been discussed in the literature). How- 
ever, also according to P&C, for samples with an adequate 
number of objects even with Eq. |] one sho uld rec over, after 



some initial fluctuations (Sylos Labini et al. |1996| ), the "cor- 
rect" dependence on distance of Eq. || for a fractal distribution, 
that is N(<R) = /<f dr 8(r -r{] ~ R D . Therefore P&C con- 
sidered an expression equivalent to Eq. [j], the integral density 
iV(< R) /R? « R D ~ 3 , and derived the value D ~ 2 from prelim- 
inary ESP data up to depths of several hundreds of Mpc. 



3. ESP survey 

The ESP final sample of galaxies (Vettolani et al. 1997, 1998), 
limited in apparent magnitude in the bj band (bj < 19.4), is 
highly complete in the fraction of measured redshifts (~ 85%). 
Despite this fact, its analysis in terms of volume-limited sub- 
samples is not entirely trivial. 
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Fig. 1. Limiting distances for ESP volume limited samples with 
different assumptions for metric and k-correction. 



N(<R) 



(1) 



where r, are the radial distances of the objects in the sample 
and 8 is the Dirac distribution. 

Now, for a simple fractal of dimension D one has 



(N(<R)) = const R D 



(2) 



One should notice an important difference between the two 
above expressions: while Eq. |l| refers to a single point (the ori- 
gin), Eq. ^ refers to a statistical property of the sample, i.e. is 



The results are somewhat dependent on the assumptions on 
the cosmological parameters, and some uncertainty is induced 
by the statistical error on the observed magnitudes (r.m.s. of 
~ 0.2 mags) and the uncertainties in the adopted k-corrections. 

In our analysis we will use the standard formula to derive 
the absolute magnitude: 



M bj 



bj-5- {5 +log[r lum (z)] +log(h)}-k(z) (3) 



where k is the k-correction term (discussed below), r\ um is the 
luminosity distance in Mpc (we assume Ho = 100 h km s 
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Fig. 2. Cumulative depth distribution of the number of ESP 
galaxies for three volume limited subsamples of case (i), cos- 
mological distance and average k-correction. The value of 
slope fitted for R > 300 Mpc/h is shown close to the relative 
line. The total number of galaxies in each subsample is given 
in the legend. 



Fig. 3. Cumulative depth distribution of the number of ESP 
galaxies for three volume limited subsamples of case (v), Eu- 
clidean distance and zero k-correction. The value of slope fit- 
ted for R > 300 Mpc/h is shown close to the relative line. The 
total number of galaxies in each subsample is given in the leg- 
end. 



Mpc^ 1 and Q.q= 1 , A = 0); henceforth we will drop the j sub- 
script from the magnitudes. 

We will consider ESP volume limited subsamples obtained 
with three cuts in absolute luminosity, namely M\, Km = —19.5, 
—20.0, —20.5. These limits span increasingly deeper volumes 
with decreasing statistics. Obviously the limiting distances of 
the volume limited subsamples are function of the adopted cos- 
mological model and k-correction. In the following we will 
consider comoving radial distance: r com — ri um / (1 + z). In the 
Euclidean case we will have r com = r\ um = cz/Hq. 

In Fig. [j] we show the limiting distances for volume-limited 
subsamples as a function of absolute magnitude for five differ- 
ent cases: 

(i) cosm ologic al distance and average k-correction as in Zucca 
et al. ( |1997| ); 

(ii) cosmological distance and k-correction estimat ed fro m the 
spectrum of each galaxy (Fiorani & Scaramella 1998 ); 

(iii) Euclidean distance and k-correction as in case (i); 

(iv) cosmological distance and zero k-correction; 

(v) Euclidean distance and zero k-correction. 

Except for case (ii) the adopted k-correction at any given red- 
shift is the same for all galaxies and in these cases the definition 
of a volume limited subsample is straightforward. In case (ii) 
we have to limit the depth in such a way to include all morpho- 
logical types, and therefore we define the limiting distance by 



using kn m (z) ~ 4.1 z, i.e. the k-correction appropriate for ellip- 
tical galaxies in our redshift range, z Si 0.25. 

For the ESP sample P&C claim that the signature of D = 
1.9 ±0.2 is seen for volume limited sub-samples at a depth 
greater than 300 Mpc/h, which according to Sylos Labini et al. 



(1996 ) is the minimum depth of statistical validity of the radial 



analysis for this sample. 

We show in Fig. || the results of fitting D in Eq. [| to the 
three volume limited subsamples of case (i) in the range R > 
300 Mpc/h. Th e influ ence of the large inhomogeneity reported 
in Zucca et al. ( 1997 ) is evident in all subsamples up to a depth 
of R ~ 300 Mpc/h, but after that we find no evidence for a slope 
~ 2, as claimed by P&C. 

Since the data presented in Fig. [| are cumulative distribu- 
tions, the points are not independent. Therefore in order to have 
an estimate of the error associated to the fit ted va lue we applied 
the bootstrap method (Efron & Tbishirani 1993 ), which yields 
a measure of the internal uncertainty of the sample at hand. 
The bootstrap estimates from 10,000 resamplings have a Gaus- 
sian shape and yield: D_i 9 5 = 2.93 ±0.14 , D_ 2 o = 3.08±0.18 
and Z)_20.5 = 3.47 ±0.28. These values are all consistent with 
D = 3, while for all of them D = 2 is at more than five standard 
deviations. If we analyze case (ii) (cosmological distance and 
k-correction estimated from the spectrum of each galaxy) we 
obtain values in good agreement with case (i), as reported in 
Table 1 . If we modify case (i) with the use of Euclidean dis- 
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tances rather than the comoving ones, we still obtain values 
which bracket D = 3 (case iii). 

Table 1 Values of D for volume-limited ESP subsamples 
Case Model k-correction M_i9.5 M-20.0 M_20.5 



1000 p 



(i) 


FRW 


(k) 


2.93 


3.08 


3.47 


(ii) 


FRW 




2.79 


3.06 


3.23 


(iii) 


Euclid 


<*> 


2.96 


2.83 


3.17 


(iv) 


FRW 





2.37 


2.37 


2.11 


(v) 


Euclid 





2.25 


1.99 


1.68 



On the contrary, if we neglect the use of k-corrections, we 
obtain values closer to D = 2 (case iv), or even in agreement 
with D — 2 (case v: D 19 . 5 = 2.25 ± 0. 10 , D_ 2 o = 1.99 ± 0.13 
and D-20.5 = 1 .68 ± 0.26), thus reproducing the result of P&C. 
The results for the last case are shown in Fig. [|. Since the effect 
of neglecting the k-correction term in Eq. |3| is to assign a lower 
intrinsic luminosity to each galaxy, the brightest subsample has 
only 43 galaxies and is dominated by discreteness. 

The results of the present analysis show the capital impor- 
tance of the k-correction term: unless one is willing to disre- 
gard the physical effect of the redshift on the observed spectrum 
and intrinsic luminosity (and moreover to adopt a purely Eu- 
clidean metric), our results do not support the claims of P&C 
of a value D = 2 for the ESP survey, but rather suggest the 
value D ~ 3. It is worthwhile to stress that the k-correction 
is an empirical and necessary correction to the magnitude of 
a galaxy, and it does not depend on any interpretation of the 
redshift. This correction simply takes into account that at dif- 
ferent redshifts we are looking at different regions of the galaxy 
spectrum. 

4. Abell clusters 

In agreement with several analyses which appeared on the sub- 



ject (s ee e.g . Bahcall 1988 an d refe rence therein), Scaramella 
et al. ( |1990| ) and Zucca et al. ( |l993| ) found that the ACO sam- 
ple suffers from a significant radial incompleteness beyond a 
depth of 300 Mpc/h. Even though the different angular selec- 
tion functions and overall completeness are different for the 
Northern (original Abell) and th e Southern sample, it has been 
argued (Scaramella et al. 199C[ ) that the radial incompleteness 
is not very strong within 300 Mpc/h. In this case, quite differ- 
ently from clustering analyses, the angular selection function 
should impact only on the amplitude but not on the radial be- 
haviour of the scalar quantity N{< R). We will therefore con- 
sider ACO clusters which have \b\ > 30" and R < 300 Mpc/h 
without any correction for angular incompleteness. 

We show in Fig. |] the results of applying Eq. || to the com- 
bined North+South sample. Not all the clusters have measured 
redshift: for those for which a measure of z is not available, 
we have used estimated z from relations which have < 20% of 
uncertainty, as reported in the papers above. We have 268 mea- 
sured z out of 363 for > 1 and 481 out of 723 for 3J. > 0. 
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Fig. 4. Cumulative depth distribution of the number of ACO 
clusters for volume limited subsamples. 



The clusters with estimated redshifts concentrate into the R 
250 Mpc/h region (within which there is ^85% of measured 
z) and can affect only the last points of Figure Q 

Repeating the same fitting as done in the ESP case, and es- 
timating similarly the errors on D, we obtain D = 3.07 ± 0.18 
and D = 2.93 ± 0.15 for richness class 11 > 1 and 1{_ > 0, re- 
spectively. 

We draw the conclusion that our analysis excludes D = 
2 with high levels of significance also for clusters, while it 
favours the canonical value D = 3. 

5. Conclusions 

On the basis of the present analysis we do not find any ev- 
idence in the ESP and ACO samples for a fractal exponent 
D ss 2 as claimed by P&C. For the ESP sample we showed 
that D ss 2 can be obtained only by neglecting the k-correction 
term: in our opinion to neglect this term would be both unphys- 
ical and unjustified. Moreover, also the cluster sample, which is 
not significantly affected by this particular uncertainty, clearly 
excludes the value D — 2. 

On the contrary, we find evidence that within the errors both 
samples show a behaviour on large scales both consistent with 
and supporting the canonical value, namely D = 3. It must be 
stressed that our results use direct depth information differently 
from indirect arguments such as scaling of angular correlation 
functions or fluctuations of cosmic backgrounds (see e.g. Pee- 



bles 1993) 



The above conclusions suggest the need of further and more 
careful and sophisticated studies on the claims of a strong ev- 
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idence for a fractal matter distribution from the data sets we 
analyzed. Also further independent checks should be done and 
possibly on other, better suited samples. 
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